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Abstract 

A mathematical theory and an accompanying numerical scheme have been 
developed for predicting the oxidation behavior of C/SiC composite structures. The 
theory is derived from the mechanics of the flow of ideal gases through a porous solid. 
Within the mathematical formulation, two diffusion mechanisms are possible; (1) the 
relative diffusion of one species with respect to the mixture, which is concentration 
gradient driven and (2) the diffusion associated with the average velocity of the gas 
mixture, which is total gas pressure gradient driven. The result of the theoretical 
formulation is a set of two coupled nonlinear differential equations written in terms of the 
oxidant and oxide partial pressures. The differential equations must be solved 
simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a 
function of space and time. The local rate of carbon oxidation is determined as a function 
of space and time using the map of the local oxidant partial vapor pressure along with the 
Arrhenius rate equation. The nonlinear differential equations are cast into matrix 
equations by applying the Bubnov-Galerkin weighted residual method, allowing for the 
solution of the differential equations numerically. The end result is a numerical scheme 
capable of determining the variation of the local carbon oxidation rates as a function of 
space and time for any arbitrary C/SiC composite structure. 
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Introduction 


The ability of carbon fiber-reinforced silicon carbide composites (C/SiC) to 
maintain its strength and stiffness at high temperatures as well as its low density make it 
an attractive candidate for many applications in future launch vehicles. These 
applications include turbomachinery components and thrust chambers in future 
propulsion systems as well as control surfaces, leading edges and thermal protection 
systems for vehicle airframes. One of the more formidable obstacles to the widespread 
use of C/SiC structures in future launch vehicles is that the carbon fibers oxidize at 
medium to high temperatures in an environment in which oxygen is present. This does 
not forbid the use of C/SiC in future launch vehicle applications, as long as it can be 
verified through testing and analysis that the component will maintain its strength and 
stiffness throughout its service life, with the demonstration of sufficient safety factors. As 
such, an assessment of the oxidation behavior of C/SiC composite structures must be 
included along with the usual design analysis activities such as the thermal, dynamic and 
thermostructural analysis of the component. It is therefore necessary to develop a tool 
that is capable of determining the spatial distribution of the extent of oxidation and the 
residual strength and stiffness in the C/SiC component as a function of the time, 
temperature and environmental oxygen concentrations to which the C/SiC structure is 
exposed. Currently, no such oxidation analysis tool is available to designers, who wish to 
utilize C/SiC composites. 

Oxygen attacks the carbon in C/SiC composites both on the surface and in the 
interior of the composite. The oxygen achieves access to the interior of the composite via 
an interconnected pore network (Lamouroux, et ah, 1993), which is formed by the 
combination of matrix cracks, separation of fibers from the pyrocarbon coating and void 
spaces between adjacent plies. The matrix cracks and fiber/coating separations are due to 
tensile stresses in the matrix, which are a result of the thermal expansion mismatch 
between the carbon fibers and the silicon carbide matrix in concert with the temperature 
excursions during processing and cool down. The large free spaces between plies are due 
to insufficient void filling during matrix infiltration. The porous nature of C/SiC 
composites is evident in fig. 1 where the microstructure of a 2-D C/SiC composite is 
shown. The oxidation of carbon in the interior of C/SiC composites is strictly tied to the 
transport of oxygen into and the transport of oxides out of the material. Any viable 
oxidation model for C/SiC composites must include the solution of species transport 
equations as the transport has a direct impact on the rate of carbon oxidation. 

Oxidation models have been developed in the past in order to study the physics of 
the oxidation process in carbon fiber-reinforced composites. Medford (1975) proposed 
one of the earliest models. His approach was used to predict the oxidation performance of 
the Space Shuttle’s carbon-carbon wing leading edge. His model attempted to simulate 
the diffusion of oxygen to the carbon-carbon substrate down a fissure in the SiC coating. 
He assumed steady state diffusion that was driven by oxygen concentration gradients 
across the coating thickness. 
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Eckel, et al. (1995) proposed a similar model to determine the oxidation reeession 
rate of a carbon fiber embedded in a non-reactive matrix. The oxygen eoncentration 
distribution in a tubular pore was determined by the solution of a steady-state gas 
diffusion equation, whieh was derived assuming that the specie migration is a relative 
diffusion in a stationary gas mixture. The oxygen eoneentration adjacent to the carbon 
fiber surface and the temperature dictate the rate of carbon surface reeession. Halbig 
(2001) adapted the model established by Eckel, et al. to simulate the fiber surfaee 
reeession in C/SiC specimens. His approaeh presupposes a eraek extending through the 
gage seetion bridged by an array of continuous carbon fibers. 

Although these studies have provided insight into the physies of earbon oxidation 
in ceramic composites, these approaches are not readily applicable to support the design 
analysis of C/SiC structures as they are impractieal for predicting the residual strength 
and stiffness as a function of space and time for any arbitrarily-shaped C/SiC structure. 
Indeed, these methods study the problem of carbon oxidation on a very fine scale, 
retaining the heterogeneous nature of C/SiC eomposites. As a result, these methods 
would be too eumbersome if they were applied to analyze C/SiC structures on a global 
level. 


The purpose of this paper is to establish and develop an aecurate analysis method 
to simulate the oxidation behavior of C/SiC composite structures in high temperature 
applications. In this paper, the mathematical theory that is the foundation of the method is 
presented. The theory is derived by assuming that the C/SiC material is a homogeneous, 
orthotropic porous body with a solid skeleton that is a mixture of multiple solid 
eonstituents where some of the constituents are reaetive. The oxygen (oxidant) and the 
oxides (produet) flow through the pore network and the partial pressures of the gases vary 
with space and time. Applying the fundamentals of porous media to this problem, namely 
the mass eonservation equation for eaeh gas specie as well as the transport meehanisms, a 
set of coupled, nonlinear differential equations is obtained. The Bubnov-Galerkin 
weighted residual method is used to cast the governing differential equations into matrix 
equations which allow us to solve the differential equations numerically. 



I .Oft •HOMI ■ I 


Eig. 1. — . Optical Microscopy of 2-D C/SiC Composite. 
(Courtesy of Southern Researeh Institute, Birmingham, AL) 
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Mathematical Theory 

Mass Continuity for Flow Through Porous Media 


In order to develop a mathematical theory to study the oxidation process, we 
make use of the basic principals of porous media theory, namely the equation for the 
continuity of mass of gaseous species flowing through a porous solid as well as 
mechanisms for transport in the porous solid. It is assumed that the oxygen and the 
products of the oxidation reaction (oxides) exist only in the pores of the material in a gas 
form and that these gases behave ideally, that the solid skeleton consists of a mixture of 
both silicon carbide and carbon in the solid form and that the carbon is oxidized at the 
interface between the solid skeleton and the pore, in other words, at the wall of the pore 
space. 


The local form of the mass continuity equations for oxygen and the oxide species, 
flowing through a porous solid body, may be written as 

+ (la) 

and 

(*^ + ?'^c„=-N2^(Mc/+3ipc). (lb) 

respectively, where (f) is the volumetric porosity, and where pP is the local partial 
density in the pore volume, M^ is the molecular weight and J. is the local mass flux 

vector for specie i. The subscripts ox, co and c refer to the oxygen, the oxide (either 
carbon monoxide or carbon dioxide) and the carbon (solid form) species, respectively. 

In equations (la) and (lb), the symbols 91 and 91 denote the local time rate 

of carbon fiber mass and pyrocarbon coating mass consumption due to the oxidation 
reaction per unit bulk volume, respectively. Further, N \ and N 2 are the stoichiometric 
constants for the oxidation reactions. That is, Nj is the ratio of the number of moles of 
oxygen consumed in the oxidation reaction to the number of moles of carbon consumed 
in the reaction and N 2 is the ratio of the number of moles of oxide produced in the 
oxidation reaction to the number of moles of carbon consumed. As such, the right hand 
side (RHS) of (la) is the local rate of oxygen mass consumed in the oxidation reaction 
per unit bulk volume and the RHS of (lb) is the local rate of oxide mass produced by the 
oxidation reaction per unit bulk volume. 
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In the oxidation of carbon, multiple reaction mechanisms are possible (Walker, et 
ah, 1959). In this study, we will concentrate our attention on two mechanisms, which are 
denoted by the balance equations as 

C + O 2 ^ CO 2 (Reaction A) 

2C + 02^ 2CO (Reaction B). 

Thus, for Reaction A, Nj = N 2 = 1 , and for Reaction B, =1/2 and N 2 = 1 . 

Upon substituting the ideal gas law for each specie i {pf = PiMj / RT ) into the 

first term in equations (la) and (lb) and upon performing the differentiation, we obtain, 
after rearranging. 


and 


^ RT dt - ^ ^ T dt 


RT dt M, 


^cf 


Pco dT 
T dt 


= 0 


(2a) 

(2b) 


Mass Flux Constitutive Relations 


The mass flux of oxygen and the oxides through the pore network occurs via two 
mechanisms: gas pressure gradient-driven flow and concentration gradient-driven flow. 
That is, we may write the mass flux vector of specie i as the sum of two mass flux vectors 
as 


J.=J.+J’’.^^ (3) 

~^i ~^i ^ ' 

where is the mass flux associated with the average velocity of the gas mixture and is 

gas pressure gradient-driven and where is the flux of specie i relative to the mixture 
average velocity and is concentration gradient-driven. 

The expression for the gas pressure gradient-driven flow is attributed to Darcy 
(1856) and is written as 

( 4 ) 

Pg 
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where k is the seeond-order material permeability tensor, is the viseosity of the gas 
mixture and p is the total gas pressure. 

The eoneentration gradient-driven flow is given by modifying Fiek’s law (Bird, 
Stewart, and Lightfoot, 1 960) as 

if (5) 

P” 

where p^is the loeal gas mixture density in the pore and is the diffusivity of gas 
speeie A with respeet to gas speeie B, and q>^ is the seeond-order areal porosity tensor. 

Making use of the ideal gas law for eaeh speeie i ( pf = PjMj / RT ) as well as the 
ideal gas law for the mixture ( = pM ^ / RT ), we ean rewrite equation (5) as 

yd ^ • Y(^) , (6) 

' Mg- p 

in order to write (5) in terms of the partial pressure pi . In (6), Mg is the moleeular 
weight of the gas mixture. 

Substituting equations (4) and (6) into (3) leads to 
io. = Dab 

p„ K1 p ^ 

(7a) 

-pL—k-Yj’co+DAB^^V-* ■YJ’co 

Pg RT p - 

and 

Ico = -pFo — k • St’co - Dab ■ Y.PCO 

p„ RT p - 

(7b) 

- /^/c, — k • V/’ox + 0.48 — <P . 

Pg RT p - 

In the derivation of equations (7a) and (7b), it was neeessary to employ Dalton’s law for 
the gas mixture, {p = Pg^ + Pco )’ the distributive property of the gradient operator 

( Yp ~ Yj’ox + YPco ) the ideal gas law for the gas mixture. 
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It is apparent that, upon substituting equations (7a) and (7b) into (2a) and (2b), we 
will arrive at two eoupled nonlinear differential equations written in terms of the partial 
pressures. These equations will be eoupled, sinee terms involving the gradients Yj’ox 
Yj’co appear in both equations. These equations will be nonlinear, sinee (7a) and 
(7b) eontain terms that involve the produet of the partial densities and the gradients of the 
partial pressures as well as terms that involve the produet of the partial pressures and the 
gradients of the partial pressures. The two nonlinear differential equations must be solved 
simultaneously at eaeh time step to obtain the partial pressures at eaeh spatial loeation. In 
the following seetions, we present the numerieal approaeh to aehieve the simultaneous 
solution of these equations. 


Determination of Oxidation Reaction Rates 


The time rate of the earbon oxidation reaetion is a funetion of the absolute 
temperature and the vapor pressure of the oxidant (Gulbransen, et ah, 1963). The 
dependenee of the reaetion rate on temperature and pressure is usually modeled using the 
Arrhenius rate equation of the form 


A 

dt 


f ^ 


V 


m 


m 


c_ 

o 

^ ''oxidation 


= -kg exp 


-Eg 

RT 


Pox 


( 8 ) 


where and are the instantaneous earbon mass and the initial earbon mass, 
respeetively, and where kg is the pre-exponential eoeffieient. Eg is the aetivation energy 
of the oxidation reaetion and n is the order of the reaetion. The Arrhenius eonstants kg , 
Eg and n are obtained by eurve fitting thermogravimetrie analysis (TGA) measurements 
to equation (8). Thermogravimetrie analysis is performed separately on the individual 
earbon eonstituents (earbon fibers and the earbon eoating). 

As we are eoneemed with a solid mixture eontaining earbon, then for any unit 
volume of material, we ean replaee the mass ratio in (8) with a density ratio and equation 
(8) ean be written as 


V ^oxidation 


= -PcK^^v 


-Eg 

RT 


Pox 


( 9 ) 


where Pg and p° are the instantaneous mass density and initial mass density of earbon 
in the mixture, respeetively. 

In applying equation (9), it is neeessary to reeognize a few key differenees 
between the oxidation of earbon fibers or pyroearbon eoating in TGA tests and the 
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oxidation of these eonstituents in C/SiC eomposites. First, the ambient partial pressure of 
oxygen in the TGA experiments is eonstant and speeified as a test eondition, whereas the 
oxygen partial pressure in the C/SiC eomposite varies with time and spatial loeation. 
Seeond, the earbon surfaee area that is exposed to oxygen vapor may be approximated as 
elose to 100 pereent in the TGA experiments, whereas this surfaee area fraetion in the 
C/SiC eomposite is eonsiderably less than this estimate and also varies as a funetion of 
time and spatial loeation. Furthermore, we note that by definition, 

(AT >’p‘ 

mass of earbon fiber and mass of pyroearbon eoating per unit bulk volume of material, 
respeetively. Given these eonsiderations, equation (9) may be rewritten separately for 
both the earbon fibers and pyroearbon eoating, as 


and 


^•«c/ =-p:f 


Pox 

\pIx) 


Ycf 


w 


cf ucf 

IV rt 


exp 


- 

RT 


(10a) 


^ pc P pc 


f 

Pox 

\pIxJ 


xifP^kP^ exp 


RT 


(10b) 


respeetively, where , Ef and are the Arrhenius eonstants assoeiated with the 

earbon fiber oxidation and kP^ , EP^ and are the Arrhenius eonstants assoeiated 

with the pyroearbon eoating oxidation. The quantity is the ambient oxygen vapor 
pressure in the TGA experiments in whieh the values of the Arrhenius eonstants were 
determined. 

Furthermore, in equations (10), we have introdueed and y/P^ as the fraetion 
of the fiber surfaee area and fraetion of pyroearbon eoating surfaee area that is exposed to 

oxygen in the eomposite material. It should be noted that y/‘^^ and y/P^ will be a 
funetion of time and the spatial loeation as it will depend on the temperature, the stress 
state and the degree to whieh the earbon is oxidized. We will leave the determination of 
the funetional dependenee of the surfaee area fraetions on the temperature, the stress state 
and the degree of oxidation as the subjeet of future studies. 

Assuming the density of earbon fiber and pyroearbon eoating are approximately 
the same, the relation between the densities and the volume fraetions are p^j- = 'pY^f 

and Pp^, = PYpc , where p^ is the density of pure earbon, is the volume fraetion of 
earbon fiber and Vp^ is the volume fraetion of pyroearbon eoating in the eomposite. 
Therefore, equation (10a) and (10b) may be written 
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Pc^cf 


( ycf 

Pox 

pIx) 


f ~cf~^ 




xj/ ■' exp 


-E 


V / 


and 


91 

pc 


Pc^ pc 


f \^pc 
Pox 

\pIx) 


xifP^kP^ exp 


'-_El_'- 

V y 


(11b) 


where and v® are the initial values of the volume fraetions. 
cj pc 


Finite Element Formulation 

We now seek a method to solve equations (2a) and (2b) numerieally to obtain the 
spatial distribution of the partial pressures for any arbitrary body at diserete times over 
any arbitrary time domain. We employ the Bubnov-Galerkin finite element method to 
oast the differential equations into a matrix form in order to perform the solution 
numerieally. In addition, we note that the nonlinear differential equations must be 
linearized within eaoh time step and solved over the time domain in a pieoewise linear 
manner. 

Applying the Bubnov-Galerkin method (Bathe, 1982), the finite element form of 
(2a) and (2b) are 


J* Ni {LHS _of _ equ.{2a)]dD^ = 0 , 

(12a) 

and 

j Ni {LHS _of_ equ.{2b)]dD^ = 0 

(12b) 

where is the domain of eaoh element and where A,- are the element shape funotions. 

For the numerioal solution of equations (2a) and (2b), we 
Pox, Pco^ ^cf 91 with the approximations Pco, 91 c/ 

replaoe the variables 
and 91 , given by 

Pox k, 0 = Ni {x, y , z)p^^^ (0 

(13a) 

Pco k, 0 = Ni (x, y , z)p^^. (0 

(13b) 

31c/ k, 0 = Ni (x, y, z)91cy. (0 

(13c) 


and 
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‘'^pc 0 = Ni (x, y, z)91 p^. (t) 


(13d) 


where p^j^- and are the element nodal values for the oxygen vapor pressure and 

oxide vapor pressure, respectively, and where 91 and 91 are the element nodal 

values for the rate of carbon fiber mass and pyrocarbon mass consumption due to 
oxidation per unit bulk volume, respectively. Therefore, using equation (12a), (12b), and 
(13a) through (13d) and noting that the shape functions are not a function of the time 
variable, the finite element form of equations (2a) and (2b) become 


I* Ni(!)^^N dD^ 
i RT J dt J 

(14a) 

- f iV,Ni^^iV,(91c/-. f =0 

J ' i yy 7 \ PCjr i p 

c M ^Pco i „ e 

[ ^-dD^ + \NiV J dD^ 

i RT J dt J ' — 

(14b) 


+ J A',N2^A'y(%^. +'y^,c,W- \ 


Using the Product Rule of differentiation along with Gauss’ Theorem (Burnett, 
1987), the second term in equations (14a) and (14b) can be expanded as 

|*A^,V-/ dD^ = iNiJ -ndY^ - \vNrJ dD^ 

J ^ox J ^ ^OX — J ^ ^ ^OX 

dC pe p,e 

and 

|*A^,V-/ dD^ = ^NiJ -ndr^- fvA^,- •/ dD^ 

J ' ^ ^CO J ' ^CO - J ^ ' ^CO 

dC pe jjC 

where T^ is the boundary of the element and, as indicated, the integration is performed 
in a closed path around the element boundary. The boundary integral terms in equations 
(15a) and (15b) allow for the application of unconstrained boundary conditions (Burnett, 


(15a) 

(15b) 


NAS A/TM— 2003 -2 12720 


10 



1987). These integral terms are only nonzero for elements where mass flux boundary 
eonditions are imposed. 

Substituting equation (7a) and (7b), equations (15a) and (15b) are written 


\NiY.-loJD^= ^Nil^^-ndr^+ \ZNi-pP^—^-Yj’oxdD^ 

dC pe jye 


+ j 


M 


ox 


RT 


Z)‘ 


^ Pco^ 

. p ) 


<p 


A 


■Yj’oxdD^+ f Wi-P^x—^-Yj^codD^ 

. Pg 


j w,- • 


Mnr 

D ab 

RT 


Z)‘ 


^ Pox^ 

VP) 


Z)‘ 




(16a) 


and 


Jw,-pP— 

d6 pe j^e 


+ I 


M 


CO 


^ Pox^ 


RT 


Z)‘ 


P ) 


(p^ • Vj^,odD^ + f VNi • pPo — k • Yj>oxdD^ 
. Pg 


(16b) 




j YNi-D - 


^AB 


RT 


Z)‘ 


^ Pco^ 

. p J 


^^■YPoxdDY 


In order to linearize equations (16a) and (16b), we treat the partial pressures 
which are operated on by the gradient operator as the solution variables and assume that 
the partial densities and the partial pressures in brackets in (16a) and (16b) are constant 
within the element over the duration of any arbitrary time step. Therefore, the partial 
pressures in brackets and the partial densities may be taken outside the volume integral 
and the gradient terms YPox YPco replaced with the approximations VA^,- 

and ^NiPco^ ■ 


Furthermore, we assume that the temperature, the gas viscosity, the diffusivity 
and the volumetric porosity are all constant within the element (although these quantities 
may vary from element to element) and therefore these quantities may be taken outside 
the volume integrals in equations (16a) and (16b). 
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After substituting equations (16a) and (16b) into (14a) and (14b), and reeognizing 
that the element nodal values of the partial pressures, and , as well as the 

temporal derivatives, dp^^. / dt and ^ functions of the spatial variables, 

equations (14a) and (14b) may be written in the matrix form 


COX 

0 

d 

o 
' — 

+ 

1 




0 

cco 

dt 

Pco j L 


KOX 

KCX 


KXC']\ Poxj 
KCO^Pcoj 


where 


[FOX'] 

\fco] 


(17) 


(COX\ \NiNjdD<^ 


{CCOlj=^^ jN;NjdD<^ 


{KOX\j 


4 ^ j Wi • H + Dab j Wi ■ ■ YNjdD^ 

/^p- RT p 


(KXC\j 


jvjv,- ■ k ■ YNjdD‘ - D^b % fvjv,- -(p -YNjdD‘ 
P CT RT p 


[KCXV = ^ f VA,- • k • VNjdD^ - Dab f VNi ■ fp 

■’ jUg \ ^ ^ RT p ^ ^ 


• VNjdD^ 


Z)‘ 


Z)‘ 


'' P 

{KCO\j = ^ j VA,- • k • VNjdD^ + Dab 


^co Pox 


Ps^ " RT p 


jvA,r 


• VNjdD^ 


Z)‘ 


(FOX), =Ni 


fx,X,(9!rf, +«„ , W)® +«»A7— {N:dD‘ 

Mi'Dcfj PCjr ^ T dt i ' 


- iXjJ -ndr^ 

J ' —ox -* 
re 


(18a) 

(18b) 

(18c) 

(18d) 

(18e) 

(181) 

(18g) 
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( 18 h) 


(FCO\ 


-N2 




- iNiJ -ndr^ 

J ' -^co - 

ye 


and where the symbol (^) above the quantity denotes the element average value of that 
quantity. 

We ean represent equation (17) more eoneisely as 

[C\jJlP\n+WA={Fl ( 19 ) 


where {p} represents the solution veetor > Pcoi J where the subseript n 

indieates the matriees and veetors at time . Using the Baekwards Differenee Method 
(Burnett, 1987), we have for the time rate of ehange of the partial pressure veetor 

4 p ]„ kL-kL-i ,20) 

dt 


where Substituting (20), we may rewrite (19) in the simple matrix form 

(21) 

where 


and 



1 

A(„ 


kL+hL 


(22a) 



1 

Atn 


kLW„-i+kL 


(22b) 


Thus, the partial pressures and p^g at eaeh node are determined by the 
solution of equation (21) at eaeh time step. The effeetive stiffness matrix J 
effeetive foree veetor for time are ealeulated using equations (18), (22a) and 

(22b) and the partial pressure values determined at the previous time step, t^_i . A flow 
ehart illustrating the steps involved in the numerieal solution routine is shown in fig. 2. 
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Fig. 2. — Flow Chart Showing the Steps Involved in the Numerieal Solution Routine. 


Discussion 

In this study, a mathematieal theory and an aoeompanying numerieal seheme have 
been developed for predieting the oxidation behavior of C/SiC eomposite struetures. The 
theory is derived from the meehanies of the flow of ideal gases through a porous solid. 
Within the mathematieal formulation, two diffusion meehanisms are possible; (1) the 
relative diffusion of one speeies with respeet to the mixture, whieh is eoneentration 
gradient driven and (2) the diffusion assoeiated with the average veloeity of the gas 
mixture, which is total gas pressure gradient driven. The result of the theoretical 
formulation is a set of two coupled nonlinear differential equations written in terms of the 
oxidant and oxide partial pressures. The differential equations must be solved 
simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a 
function of space and time. The local rate of carbon oxidation is determined as a function 
of space and time using the map of the local oxidant partial vapor pressure along with the 
Arrhenius rate equation. The nonlinear differential equations are cast into matrix 
equations by applying the Bubnov-Galerkin weighted residual method, allowing for the 
solution of the differential equations numerically. The Backward Difference Method has 
been employed as the time marching scheme. As the differential equations are nonlinear 
they must be linearized within the finite element formulation. The nonlinear differential 
equations are linearized within each time step and solved over the time domain in a 
piecewise linear manner. This is achieved by continuously updating the system stiffness 
matrix and system force vector based on the values of the solution variables determined 
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in the previous time step. The end result is a numerieal seheme eapable of determining 
the variation of the loeal earbon oxidation rates as a funetion of spaee and time for any 
arbitrary C/SiC eomposite strueture. 

The differenee between the modeling approaeh proposed here and those 
mentioned previously is that the eurrent method uses a eontinuum meehanies approaeh to 
modeling the oxidation proeess in C/SiC eomposites. This allows for the determination of 
the earbon oxidation rate and earbon volume fraetion distributions on a global level for 
any arbitrarily-shaped C/SiC eomponent. Furthermore, the eurrent approaeh ineludes an 
additional diffusion meehanism whieh is driven by total gas pressure gradients. This 
approaeh will prediet a signifieantly different diffusion behavior (and eonsequently a 
different oxidation pattern) than an approaeh that negleets the Darey flow meehanism 
whenever the total pressure gradients beeome signifieant. One instanee where total 
pressure gradients may be signifieant is in C/SiC struetures where high thermal gradients 
are present and during transient heating eonditions. Finally, the eurrent approach allows 
for a transient solution to the diffusion equation, whereas the previous methods used only 
a steady-state diffusion solution. Obviously, the current approach is a significant 
improvement in this regard, as the diffusion conditions present in an oxidizing C/SiC 
composite are almost always transient. 
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